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Abstract 



o 
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^ ' In this work we give a mechanical (Hamiltonian) interpretation of the so called 

■ spectrality property introduced by Sklyanin and Kuznetsov in the context of Back- 

lund transformations (BTs) for finite dimensional integrable systems. The prop- 
erty turns out to be deeply connected with the Hamilton-Jacobi separation of 
' variables and can lead to the explicit integration of the underlying model through 

, the expression of the BTs. Once such construction is given, it is shown, in a sim- 

ple example, that it is possible to interpret the Baxter Q operator defining the 
quantum BTs us the Green's function, or propagator, of the time dependent 
Schrodinger equation for the interpolating Hamiltonian. 
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^ I maps. Quantum propagator 

1 Introduction 

^ ' Starting from the last decade of the previous millennium, a number of results on 

■ discretization of finite dimensional integrable systems appeared: we quote for instance 

the Euler top [1], the Lagrange top [2], the rational Gaudin magnet [5] the Ruijsenaars- 
Schneider model [10], the Henon-Heiles, Garnier and Neumann systems [12], [13], [4] 
and others (see also [17] and the references therein). It turned out [8] that all the 
exact discretizations for these systems associate new solutions to a given one: they are 
BTs for such systems. These developments suggested to Sklyanin and Kuznetsov [6] 
that the concept of Backlund transformations should be revised in order to enlighten 
some new (and old) aspects of the subject. Actually, these two authors had previously 
obtained fundamental results about "separation of variables" and their connection 
with the new techniques of classical and quantum inverse scattering method. These 
findings paved the way to a better understanding of the role of BTs in the context 



^ "This paper is a contribution to Solitons in 1+1 and 2+1 dimensions. DS, KP and all that, 
conference in honor of the 70th birthday of Marco Boiti and Flora Pcmpinelli, Lccce, 2011 September 
13-14." 
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of finite dimensional integrable systems and to unveiling their deep and fruitful links 
with Hamiltonian dynamics and separation of variables (for a characterization of the 
BTs for finite dimensional integrable systems by a geometrical point of view see [8]). 
Let us briefly remember some of these findings. Assume that the dynamical system 
under consideration is defined by a Lax pair with a spectral parameter, say L{X) 
and M(A). The trace of the powers of L{X) will be conserved quantities, so also 
the determinant of the Lax matrix will be conserved. The BTs, in the case of finite 
dimensional systems, are canonical transformations preserving the algebraic form of 
the integrals, so if {pi,qi)^=i is the set of the dynamical variables, one has: 

(Pi,9i)r=i =^ (Pi'^OLi L{X,pk,qk) ^ L{X) = L{X,pk,qk) (1) 

{Pi,Pj} = {Pi,Pj}; {Qi, Qj} = {Qi, Qj}; {Pi, Qj} = {Pi, Qj} = ^ij- (2) 

Since the algebraic form of the integrals is preserved, also the determinant of the Lax 
matrices L{X) and L{X) is the same: 

Hi{pk,qk) = Hi{pk,qk) =^ det(L(A)) = det(L(A)) 

Hence, the existence of a BTs for the given system entails that the two Lax matrices 
L(A) and L{X) are connected by a similarity transformation, provided by the so called 
dressing matrix -D(A): 

L{X)D{X) = D{X)L{X) (3) 

Obviously the dressing matrix is not unique, because one can have different BTs for the 
same system. Furthermore, if the transformations are not explicit, -D(A) might depend 
on both the sets and {Pi,qi)^=i More interesting are parametric and explicit 

BTs, that is BTs depending on one or more parameters, arbitrary (or restricted to some 
range of values) such that the expression of the variables {pi, is known explicitly 
in terms of the variables {pi,qi)i'=i- A sufficient condition to obtain parametric and 
explicit BTs is to find a dressing matrix D{X) that i) whose manifold coincides with a 
symplectic leaf of the same Poisson bracket as those satisfied by L{X) [TB], and ii) such 

that its determinant has a non dynamical zero , i.e. det {D{X)) = , where fi is 

a parameter not depending on the dynamical variables. The condition i) ensures that 
the transformations are canonical (see [H]), the condition ii) ensures that are explicit 
and parametric. Suppose indeed that the particular symplectic leaf defined by -D(A) is 
one dimensional, so that -D(A) will depend on another (up to now free) variable, say a. 
Then, since det(D(/i)) = 0, D{fi) will have a kernel, say that in general also will 
depend on /i and a. But from [3] we see also that \Q) is an eigenvector of L(/i); indeed: 

L(/i)fz}(/i)|fi(^))) =D(^)fL(^)|fi(/i))) ^L(/i)|fi(^))=7(/i)|l^(^)) (4) 



The last equivalence can be read as an equation for a that, if solved, will give to a a 
dependence on /i and on the dynamical variables of the system (only one set of them, 
in this case the "untilded one"). So, returning to [H one obtains parametric (/x) and 
explicit transformations. 
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Depending on the dimension of the leaf defined by D{X) it is possible to obtain multi- 
parametric BTs. Also, one can obtain multi-parametric BTs through repeated itera- 
tions of one parameter transformations, since BTs for finite dimensional systems always 
commute [TB] . 

For simplicity in the rest of the paper we take L{\) to be a st(2) matrix, also if our 
results can be easily generalized to L{X) G 0t(2). For BTs found through a NxN Dar- 
boux matrix see [10], [16]. Please note that in the case L{X) G sl{2) the function 7(/i) 
in m satisfies 7^(/i) -l- det(L(/i)) = and is the generating function of the conserved 
quantities of the system. 

An additional property of BTs was introduced in ^ , namely the "spectrality property" . 
In order to understand this property remember that BTs are canonical transformations, 
implying that it exists a generating function q, /i) such that: 

Pi = -K— Pi = -77^ 

OQi OQi 

The spectrality property [6] says that there exists a function /i satisfying the following 
equation: 

= -/i(7(^),/i) (5) 

where 7(//) is just the function in|4) it satisfies the relation det(L(/i) — r;(/i)ll) = 0, that 
can be seen as the separation equation, in the sense of Hamilton- Jacobi separability, 
for the dynamical system; by a classical point of view one obtains the quadrature of 
the equations, by a quantum point of view the factorization of the eigenfunctions [ISj 
[6], [8]. 

Starting by this point we show how the spectrality property can lead to the "lineariza- 
tion" of the maps defined by the BTs, that, in turns, give the general solution of the 
equations of motion. This construction give also us the possibilty to properly interpret 
the Baxter Q operator (representing the quantum BTs) as the Green's function of the 
Schrodinger equation defined by the Hamiltonian interpolating the discrete flow given 
by the BTs. 



2 Spectrality property and separation of variables. 

First of all we want to make some comments about the equation |5] defining the spectral- 
ity property. Suppose that it is possible to choose the parameter /i in such a manner 
that the dressing matrix D{X, fi) is proportional to the identity when fi = 0. Expanding 
equation [3] around = 0, if D{X,fi) = /c(ll + fiDo{X) + 0{fi^)), one obtains: 

L(A) = L(A) + /i[Do(A), L(A)] + Oifi') (6) 

so that in the limit — ?■ 0, by defining L = lim^^o "^^5 BTs define an Hamiltonian 
flow: in this sense n can be considered an evolution parameter. Equation |5] than 
resembles the Hamilton- Jacobi equation with respect to the time since the function 
7(/.i) contains all the conserved quantities of the system. This point will be made 
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clearer in the next lines. 

Let us assume to have a set of BTs, with the parameter /u playing the role of a time, 
in the sense given before: 



Pi= Pi{Pk,qk,fJ-) with Pi 
qi = qiiPk,qk,IJ^) with qA 



tj.=0 



Pi 



(7) 



These transformations can be also rewritten as (assuming obviously the implicit func- 
tion theorem can be applied): 

Pi = Pi{qk,qk,fj') 

Pi = Pi{qk,qk,fJ') 

Since [7] and [H are canonical transformations, there exists the respective generating func- 
tions, say J^Q{p,q,fi) and solving the corresponding system of differential 
equations: 



Pi 



^ dqk dJ^o 



E~ <jqk 



k=l 



~ oqk 



k=l 



dqi dqi 
dpi 



Pi 



Pi 



dqi 
dqi 



(9) 



Now we assume that the transformations possess the spectrality property both with 
respect to J-o(p, g, ^) and with respect to g, So there exist two functions, say 
/ and (7, such that: 



5Jo 

9yU 



-/(7(/^),/^) 



(10) 



This assumption has non trivial consequences as will be shown in the following two the- 
orems. The result will be that, under the assumption, the BTs can be re-parametrized 
so to represent the solution of the Hamilton- Jacobi equation for the Hamiltonian in- 
terpolating the flow defined by [61 

Theorem 1 Suppose that the spectrality property holds true both forj^i and J-q. Then 
Y.k=iPk^ 9^'"^^ ^y^~ ^^^^ Y2=iPk^ « function of only and /i. 



Proof. The two generating functions J-q and are related by [3]: 

J^oip, q, fj-) = J^i{q{p, q, /^), g, /^) 

So, utilizing equation [TT] and pi = — it is possible to write: 



:ii) 



djj, 



k=l 



'k dfi dpi 



q=const. 



n 



k=l 



dpi 



□ 
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Theorem 2 Suppose that the assumption of Theorem [I] holds, then giving to a 
dependence on the constants of motion through any of the root(s) of the equation 
/i = /i((yf(/i, 7(/i))), where h is an arbitrary function, one obtains again a set of BTs. 
The generating function of the new transformations is given by: 



F{p,q,^i^^^) = J'o{p,q,fi 



g{-f{fi),fi)dn 



(12) 



where jj,^^"^ = h{g{jJL^^\'y{iJL^^'>))) is any of the root of fi = 7(/i))) 



Proof. Suppose to give to /x a dependence to the dynamical variables {pi, qi)i'=i- We ask 
if it is possible to choose some function q) in such a way that the transformations 
Pi = Pi{P^Qi fJ'iPiQ)) and g,j = qi{p,q, fi{p,q)) obtained inserting the function fi{p,q) in 
the expressions [8] are again canonical. 

If the new transformations are canonical they possess a generating function, say J-", 
such that: 



Pi 



k=l 



sr^ ~ dqk _ 
^^^^ dqi dqi 

- 

dpi 



\ ^ . oqk 

^ dpi 

k=l 



(13) 



Inserting the ansatz J-" = J^Q + A{p, q) into the previous system and taking into account 

dfj. 

to be satisfied by A{p, q): 



that by Theorem [T] 51 = — Yll:=iPk^^ readily finds the equations that has 



dA 


dii 


dqi 


= Sir 

dqi 


OA 




dpi 


OPi 



l..n 
l..n 



Indeed from T = 



A and from [13] we have: 



dpi dpi 



IJ,=const. 



dTo dfi dA ^ dqk 

dji dpi dpi ^^^^ dpi 



fj.=const. 



n 

Y^Pk 

k=l 



'k dfi 



dfi dpi 



(14) 



(15) 



When fi is constant the BTs are canonical transformations with generating function 
J-Q, so it holds the equation: 



dTo 



dpi 



fj,=const. 



- dqk 
dpi 



Inserting this equivalence in 



k=l 

we are left with: 



li=const. 



dJ^o , ~ dqk\ dfi dA 

2_^Pk— I — + — = 



9/i 



fc=i 



dpL j dpi dpi 
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But from Theorem [T] we have that + Ylt=iPk^ = ~9i so finally: 

dA dfi 
opi dpi 

In the same manner, by the equation for |^ it is possible to show that: 

dA dfi 

7r = 9^ i = l..n 

oqi dqi 

The compatibility equations for the function A{p,q), that is ^gA_ = give (by 

dg_di^ ^ dg_dfi_ .^Q. 
dpk dqi dqi dpk 

These equations means that /i has to be an arbitrary function of g, that is: 



Note that this (set of) implicit equation for fi fixes the dependence of /i by the dynamical 
variables {pi,qi)i'=i through its roots. At any of the root of this (set of) equation will 
correspond a function g and a function /i, both depending only on the dynamical 
variables {pi,qi)^=i, satisfying the set of differential equations [161 Indicating the k^^ 
root of /i = h{g{ii,'j{fi))) as fi^'''^ and returning to the system [HI one finds: 

MPi(1)= 9{l{f^),f^)dfi (17) 

□ 

Let us now make some remarks that will be useful in what follows. 

Remark 1 The generating function J^{p, q, fi) as a function of fi satisfies the equation: 

dTjp, q, ji) ^ _ 

dji 5/i 

This is just a trivial calculation on ^-^^^'^'^^ for the function J^{p, q, /i) as given in [12] 
and taking into account the Theorem [H 

Remark 2 By adding a dependence on an extra constant parameter, say T , to the 
function fx, that is by considering fi = ijS^\T), the canonicity of the transformations is 
preserved, so again one has a set of BTs. 

This remark gives the possibility to obtain a new parametric family of BTs by letting 
to (jS^^ to depend also by the parameter T. Indeed it is possible to repeat the same line 
of reasoning of the Theorem [2] by considering, from the beginning, as a function of 
the dynamical variables (pi, ^i)"=i and the new constant parameter T. Obviously now 
also the generating function J-" will got a dependence on the parameter T. 
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The key point of all this construction is that, assuming the spectrality property holds 
both for J-Q and J-*!, it is possible, at least in principle, to obtain a new, larger, family 
of parametric BTs by giving to /i a dependence on the constants of motion and on the 
parameter T: as we will show now, to this freedom to obtain new BTs it corresponds 
the possibility to get, with a suitable choice of the function n{p,q,T), the canonical 
transformation from the variables {pi,qi)i=i to the new variables {pi,qi)^=i such that 
the Hamilton- Jacobi equation, where the parameter T plays the role of "time" and the 
interpolating Hamiltonian the role of its conjugated variable, is identically solved. 
Summarizing, given a parametric BTs {pi{pk,qki fJ'),Pi{Pk,qk, IJ-)) having the spectral- 
ity property both for J^q and J-'i, the goal would be to find a function fi fitting the 
Hamiltonian- Jacobi equation for the flux described by the BTs: 

where with the italic T-L we mean the interpolating Hamiltonian of the BTs. By Remark 
[T]the previous equation can be better rewritten as: 

fc=i 

Indeed now by Theorem [1] the sum Y12=i ^ function of /i and of the dynamical 

variables only through the constants of motion; obviously also the interpolating Hamil- 
tonian "H will depend only on some combinations of the constant of motions. So the 
equation [18] is indeed solvable in terms of /i. Inserting the solution fi{p, q, T) into the 
expressions of BTs we obtain, by construction, the general solution of the equations of 
motion with respect the interpolating Hamiltonian "H. 

Let us make another remark. It seems that there is an ambiguity in the choice of the 
function fi solving equation [I8] because, of course, the solution depends also on an ar- 
bitrary function of the constants of motion. But now we will show that the freedom to 
choose the value of this arbitrary function corresponds to the freedom in the choice of 
the initial conditions in the general solution of our equations of motion: so it represents 
a shift in the time T. Indeed, taking into account Theorem [1] and the equations [101 
the solution of the Hamilton- Jacobi equation [TS] can be implicitly written as: 



(19) 



where, by [lO] we recall that the functions / and g are respectively given by 
and and F is an arbitrary function of the Hamiltonians of the system Hi. It is 

now evident from equation [TH] that i) the arbitrary function F appears as an additive 
quantity with respect to the time parameter T and ii) if one wants to retain the initial 
conditions as given in [TJ that is Pi\j.^Q = Pi and qi\rp^Q = qi-, then one has to choose 
F = 0. 

In the next section we will give a simple application of the theorems just seen with the 
BTs for the one dimensional harmonic oscillator. 
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3 A tutorial example: the harmonic oscillator 

A Lax representation with spectral parameter for the one dimensional harmonic oscil- 
lator is given by: 

m=(^,l, M = ^(J L{X) = [L,M] 

In this simple case the spectral curve is simply given by: 

2 2 

f(A) = -det(L(A)) = l + ^^ 

In the following for brevity we pose a = p — iq and a* = p + iq. 

A class of BTs can be obtained by a dressing matrix D{X) parametrized as follows: 



D{X) 



1 f 
f 1 



As pointed out in the introduction, to find explicit transformations we impose that 
det{D{X = fi)) = 0. This constrain means that a/3 = yU.^, so we can rewrite, posing 

1 ^ 



DiX) 



EL I 



A 

Obviously when X = fi the dressing matrix possesses a kernel Explicitly it is 

given by: 

Then the eigenvector relation: 

gives C as a function of a and a*: 

/^(l-7(^)) a* 



C 



a /i(l + 7(/i)) 

Now we are able to write out the explicit form of BTs: 

7(/i) + 1 



a = a- , , _ 2 



a = a 



^'"/-^ w,th 7M' = 1 + ^, 

7(/^) + 1 



2 



Note that when /i is equal to zero the transformations reduce to the identity map. 
What we want to check now is whether these canonical transformations possess the 
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spectrality property both for J-q and for T\. The generating function in the variables 
a, a* is given by J-'o{a, a* , fi) = —2jj^{j{fi) + 1), while the inverse transformations: 



2^* 



-,*\2 



2~* 



,*^2 



are generated by the function J^i{a* , a* , fi) = ff_", . By differentiating J-'q and J-'i one 
finds: 



, 7(m)-1 
7(m) + 1 



-4^(7(/i) H 



so that both and are only functions of the integral {p^+q^) and of the parameter 
^. The results of the previous section can be applied. In the one dimensional case the 
Theorem [2] simply implies that /i can get an arbitrary dependence on the Hamiltonian 
and the BTs are again canonical transformations. So by now we consider /i as a function 
of |a| and of a new parameter T, that is ^ = //(|a|, T). Again from Theorem |2] the new 
transformations are again canonical with the generating function: 
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^l=^l{\a\) 



In order to fix the dependence of by \a\ and T we need to find the interpolating Hamil- 
tonian. As explained at the beginning of the section [2] the interpolating Hamiltonian 
of the flow with respect to /i can be extrapolated considering the relation LD = DL 
in the limit — j- 0: 



D = l-fiDo + 0{n^) L(A) = lim 



L-L 



(20) 



In our case the matrix -Do (A) is given by: 



DoiX) 



A \ a 



The continuous flow given by|20]is governed by the Hamiltonian "H = 2i|a|. Note that 

the factor "i" appear because a and a* are not canonically conjugate. 

Now we have all the ingredients of the recipe. So we can look at the Hamilton- Jacobi 

equation: 



ri — a^:^^:^ = 



dT 

— )■ fi{\a\,T) = — i|a| sin 



2i|a| + 
T 



2 


a\ 


2 


djj 




a 


2dT 







(21) 
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As noted at the end of section [2], one can choose the arbitrary function -F(|a|) to be 
zero by requiring that fi 



= 0. 



Inserting the formula n{\a\,T) = — i|a| sin(|^) into the expressions of BTs [3] and re- 



turning to the physical variables p = ^-±2 and q = one finds: 



^ ^ 2T 2T 

P = P cos( =) - q sin( 



< 



2T 2T 
q = gcos( — ) +psin( 



that is the general solution of the equations of motion governed by the Hamiltonian "H = 
2\/p^ + q^- Note that we now omit the factor "i" from the interpolating Hamiltonian 
because we are considering directly the flow with respect to T and not with respect to 
fi (when T goes to zero /i ~ — iT). 

Note that the parameter T is a "linearizing parameter" for the flow. 
For completeness we have to cite [T^ , where another way to find the functional relation 
between fi and T (if any!) has been given. However in [H] there were no definite 
statement about the necessary and sufficient condition for such relation to exist, neither 
its connections with analytic mechanics was pointed out. Nevertheless that point 
of view can be very useful because to the explicit knowledge of the BTs does not 
always correspond the explicit knowledge of the generating functions J-q and of the 
transformations. It can be also shown that applying the machinery described in |14j . 
one obtains exactly the result /i(|a|, T) = —i\a\ sin(^), as in[2Tl The reader is referred 
to this work for other examples of linearization of BTs. 



4 Quantum Backlund transformations: some hints 
and insights. 

By a classical point of view the BTs associate new solutions of the equations of motion 
to a given one. So, in general, if g{p{t), q{t)) is a physical observable connected to the 
curve in the phase space {p(t),q{t)), g{p{t),q{t)) will be the same physical observable 
connected to the curve(p(t), g(t)). 

By a quantum point of view the BTs are represented by a unitary operator Q^^-* (// 
being the (set of) parameter(s) of the BTs) realized as an integral operator on the 
space of eigenfunctions [H], [6]: 

Q«:V'(g)^ j f^'\q,qmq)dq (22) 

The similarity transformations induced by Q^/'* are the equivalent of the classical canon- 
ical transformations and the kernel f^^\q,q) is given, in the semiclassical approxima- 
tion, by: 

/(^)(g,g)~exp(-^Fi(g,g)), h^O (23) 
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As far as we know, the explicit construction of the Q operator, known in hterature 
as Baxter operator, and its relations with BTs, have been pointed out only for the 
Toda lattice [11] , [16] and for the discrete self-trapping (DST) model [7] . By our point 
of view, when the BTs can lead to the integration of the equations through the con- 
struction given in the previous sections, the kernel f^^\q,q) has to be identified with 
the propagator or the Green's function for the time dependent Schrodinger equation, 
the Hamiltonian being that interpolating the flow represented by the BTs. Indeed, in 
this case, to the classical time shift on the trajectory specified by the initial conditions 
represented by the BTs it corresponds, as it is clear also from eq. [22l a quantum prob- 
ability amplitude for the particle with position eigenvalue q at the time to be found 
at a later time ^ in q. 

In the following we make a check on the Baxter operator for the harmonic oscillator 
by showing that indeed its kernel gives exactly the well-known propagator of the cor- 
responding Schrodinger equation. In order to compare the results we make a slight 
modification in the choice of the function fi in by posing T = Now play 

the role of the time. With this choice indeed the interpolating Hamiltonian is just the 
physical one, and obviously the BTs give the formulae: 

p = p cos(</)) -|- q sin(0) 
q = q cos(0) — p sin(0) 

In order to find the generating function Fi {q, q) we express p and p in terms of q and 

q — gcos(0) 



p 
p 



sin(0) 
q cos(0) — q 



sin(0) 

From dFi = pdq — pdq one readily finds: 

2sm(0) 

where is an arbitrary function of 0. 

Now we have to find the Q operator, or, that is the same, its kernel f^{q,q). In the 
(bra|ket) notation the kernel f^{q, q) is given, in the q representation, by {q\q). So, by 
the following form of the BTs: 

q = q cos(0) -|- p sin(0) 
q = q cos(0) — p sin(0) 



it is easily found that f {q,q) solves the system: 



gj(i) = g~j(i)cos(0)-in^sin(0) 
qf*W = qf*W cos((/>) + i^^^ sin(0) 
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where " means complex conjugation. The solution of the previous system is: 



i ( {q^ + q^) cos(0) — 2qq 



f'-'\q,q) = C{c^)exp 

\ ri \ 2 sm(0) 

in exact agreement with [231 

For completeness we recall the classical results on the propagator of the Schrodinger 
equation for the harmonic oscillator. In general the propagator can be formally written 
as: 

K{q,q) = {q\exp{ —)\q) = ^{q\n){n\q) exp( ^) 

n 

In the case of the harmonic oscillator {q\n) are the Hermite polynomials, so that: 



\ 2h J 



<1 ^TT r ^ ^e~'^"+^J 



It is possible to sum up the series by using the so called "Meheler's formula" (see e.g. 

my- 

obtaining: 

^ ' ^/MHM^ V 2sin(0) 
This formula agrees with that of f^^\q,q) identifying C(0) with , ^ =. 

A/27riSsin(0) 

On the same line of reasoning we can also ask for the quantum probability amplitude 
for the particle with momentum eigenvalue p at the time to possess at a later time fi 
the momentum p. Again there is the correspondence between the classical generating 
function, now given by dj^2 = qdp — qdp, and the kernel of the corresponding Baxter 
operator Q^^^: 

f^'\p,p) ~ exp(--^F2(p,p)), h ^ (25) 



From the equations: 



one readily finds: 



pcos(</)) — p 

sin(0) 
p — pcos(0) 
sin(0) 

. (p2 +p2)cos(0) - 2pp 
2 sm(0) 



where again r] is an arbitrary function. Whereas from the following form of the BTs: 

p = p cos(0) + q sin(0) 
p = pcos{(j)) — gsin(0) 
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one finds, recalling that in the p representation q = i^^, the equations to be satisfied 
by as: 

= p/''^'* cos(</)) — — sin(0) 
op 

/)f*(2) 

pj*{2) ^ pj*{2) pQg(^) ^ i^^^ sin((^) 
so that: 

f^'\P,p) = V{<P)e.p (4 ( (p'+P')--<4>)-2pp 
again in exact agreement with 



h \ 2sin(0) 



5 Conclusions 

In this work we tried to shed some light on the spectrality property of BTs, by showing 
how it can lead, through the constructions of the section [31 to the explicit integration 
of the underlying equations of motion. It could be interesting to apply these results 
to many-body systems possessing known BTs with spectrality property, such as, for 
example, the Toda lattice. Also the interpretation of the Baxter operator for such 
systems as the propagator for the corresponding Schrodinger operator could lead to 
interesting summation formulae for orthogonal polynomials, similar to the Meheler's 
formula [Ml As a final remark we point out that the non-obvious relation among Back- 
lund transformations and the Green's function of the Schrodinger operator provides 
a bridge between the theory of such transformations and the Feynman path integral 
approach, opening unexplored perspectives in statistical mechanics and quantum field 
theory. 
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